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Abstract — This paper addresses coordinated downlink beamform- 
ing optimization in multicell time-division duplex (TDD) systems 
^vi where a small number of parameters are exchanged between cells 
I but with no data sharing. With the goal to reach the point on 
..the Pareto boundary with max-min rate fairness, we first develop 
— >a two-step centralized optimization algorithm to design the joint 
beamforming vectors. This algorithm can achieve a further sum-rate 
^improvement over the max-min optimal performance, and is shown 
i to guarantee max-min Pareto optimality for scenarios with two base 
stations (BSs) each serving a single user. To realize a distributed 
solution with limited intercell communication, we then propose an 
iterative algorithm by exploiting an approximate uplink-downlink 
duality, in which only a small number of positive scalars are shared 
i between cells in each iteration. Simulation results show that the 
proposed distributed solution achieves a fairness rate performance 
close to the centralized algorithm while it has a better sum-rate 
C/3 performance, and demonstrates a better tradeoff between sum-rate 
Oand fairness than the Nash Bargaining solution especially at high 
1 signal-to-noise ratio. 

t-H \ Index Terms — Distributed processing, Multicell downlink beam- 
forming, Pareto boundary, Uplink-downlink duality. 

m ; 

I. Introduction 

OO . 

i , Increased frequency reuse has been an attractive solution to 
provide coverage and accommodate more users in cellular net- 
works with precious spectrum resources. In conventional cellular 
{Nl networks based on single-cell processing, the users in a cell are 
■^Honly served by the base station (BS) centered in that cell and 
|>a BS is not supposed to assist in enhancing the channel links 
'"jnot belonging to its cell, while the intercell interference (ICI) is 
dealt with by careful frequency planning and usually treated as 
j^j noise at the receiver. However, in this single-cell approach, the 
cell-edge users tend to suffer severe ICI and exhibit much poorer 
performance than the users in the interior of the cell. 

To overcome this, multicell cooperation between BSs has re- 
cently emerged as a promising solution, which has been exten- 
sively studied in both academic research and standard bodies such 
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as LTE- Advanced/IEEE 802.16m. The idea is that neighboring 
cells form a cooperating cluster by exchanging information with 
each other. To be practically useful, only a limited number 
of cooperating cells should be clustered for the coordinated 
transmission/reception. In this regard, several solutions including 
network-centric [1], user-centric methods [2^4], or combinations 
thereof [5] have been developed recently. The static network- 
centric method is simple for implementation but users at the 
cluster edges would still suffer severe out-of-cluster interference. 
This motivated the development of dynamic user-centric clustering 
methods [2^1], which allow the cell clusters to be partially 
overlapping, resulting in a significant increase in scheduling 
complexity. To tradeoff between performance and complexity, a 
combination of these two methods was recently proposed [5]. It is 
realized by using the user-centric method within a set of candidate 
clusters predetermined via a network-centric means, and is widely 
adopted by standard bodies to study the multicell cooperation. 

Another issue for the cooperative multicell system design is how 
to optimize the cooperative transmission/reception between the 
cells. Early work in the area largely focused on the uplink [6,7], 
where several BSs jointly decode the signals from the users served 
by these cells. The problem of downlink multicell processing, 
however, is much more challenging, especially with multiple 
antennas at each BS. If a joint transmission strategy is employed, 
in which both data and channel state information (CSI) are shared 
between the BSs, the cooperating cluster will essentially collapse 
into a single-cell multiuser multi-antenna system. Its potential 
gain has been investigated in [8,9]. Nevertheless, it is noted that 
per-BS power constraints should be carefully addressed and an 
optimal joint transmission requires full phase coherence between 
the signals from different BSs, which is difficult to realize in 
practice. Moreover, the required signaling between the BSs is 
vast, making it challenging for practical applications. Therefore, a 
distributed solution with limited intercell communication is highly 
desirable. 

More recently, the achievable rate performance of downlink 
multicell processing in decentralized form [10, 11] and with lim- 
ited backhaul capacity [12] was analyzed, and showed a large po- 
tential. Also, several schemes for distributed multicell processing 
have been developed in [13-15]. In particular, given the avail- 
ability of data sharing, [14] addressed the design of distributed 
transmit beamforming by recasting the downlink beamforming 
problem into a linear minimum mean-square-error (LMMSE) 
estimation problem. However, the requirement of information 
exchange between the BSs is still too high and global convergence 
is not guaranteed. Later, [15] devised a distributed design using 
only local CSI, and demonstrated that performance close to the 
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Pareto boundary can be obtained. The main issue for [14, 15] is 
that both need data sharing between the BSs, which is not likely 
to be viable due to the limited backhaul capacity. 

Motivated by this, [16-20] studied coordinated beamforming 
techniques without data sharing. In [16], Lagrangian duality theory 
was applied to design coordinated beamforming with limited 
intercell information exchange but the solution only applies to 
the sum-power minimization problem with a signal-to-interference 
plus noise ratio (SINR) constraint. Also, restricting to limited 
intercell communication, [17] presented a hierarchical iterative 
solution to maximize the worst-user rate, but global optimality 
is only obtained in the case of a total BS power constraint and 
extensions to per-BS power constraints are nontrivial. It is also 
possible to devise distributed schemes for coordinated precoding 
based on the signal-to-generated-interference plus noise ratio 
(SGINR) [18] and the signal-to-caused-interference ratio (SCIR) 
[19], with each BS performing optimization only based on local 
CSI. However, in most cases this strategy is only sum-rate optimal 
in the high signal-to-noise ratio (SNR) regime with two BSs. 

In this paper, our objective is to address the distributed real- 
ization of coordinated multicell multiuser downlink beamforming 
with limited intercell communication. Considering per-BS power 
constraints, it is known that the optimal rate region achieved by a 
coordinated multicell system is characterized by the strong Pareto 
boundary [11], where it is impossible to improve the rate of any 
user without simultaneously decreasing at least one of the others. 
With the goal to achieve a fairness based rate optimality, we target 
at the particular point on the strong Pareto boundary that has the 
maximum worst-case rate, i.e., the Pareto optimality with max-min 
fairness. Different from the conventional max-min criterion which 
is only weakly Pareto optimal 1 and usually produces identical user 
rates, this objective could achieve a good tradeoff between the 
sum-rate and user fairness. By recasting the original downlink 
optimization problem into an approximate dual virtual uplink 
problem, we will propose an iterative algorithm to reach this 
desired point. 

The main contributions of this paper are as follows: 
« We first devise a centralized solution to achieve the Pareto 
boundary with the max-min rate fairness, which is realized in 
two steps: (1) a weighted SINR balancing optimization and 
(2) a Pareto improvement. We prove that our method is able 
to reach the point of the Pareto boundary on which the worst 
rate among users is not less than any other points, especially 
for the case of 2 cooperating BSs. 
• Then, we propose a distributed algorithm to approach the 
performance of the centralized scheme, in which only a 
small number of positive scalars are exchanged between the 
BSs. The distributed solution is realized by exploiting an 
approximate uplink-downlink SINR duality subject to per- 
BS power constraints. The idea is to iteratively perform 
the weighted SINR balancing optimization in a distributed 
fashion by converting it into a virtual uplink problem. Global 
convergence of the proposed iterative algorithm is proved. 
Our simulation results also show that the proposed distributed 

'A point is called weakly Pareto optimal if it is impossible to strictly improve 
the rates of all the users. If not explicitly stated, the Pareto optimality used in this 
paper refers to strong Pareto optimality. 



solution obtains a worst-user rate performance close to the 
optimal centralized scheme while has a better sum-rate per- 
formance; and comparable sum-rate performance as the Nash 
bargaining solution (NBS) but with better max-min fairness. 

This remainder of this paper is organized as follows. Section 
II describes the system model for coordinated beamforming. 
Section III presents a two-step centralized multicell downlink 
beamforming scheme to reach the Pareto boundary with max-min 
fairness. Section IV analyzes the uplink-downlink SINR duality 
with per-BS power constraints and proposes a distributed solution. 
Simulation results are provided in Section V and we conclude the 
paper in Section VI. 



II. System Model 

Consider a multicell downlink system where cells are al- 
lowed to perform coordinated transmission. As shown in LTE- 
Advanced standardization work for coordinated multi-point trans- 
mission/reception (CoMP), it is efficient to dynamically determine 
the coordinating cell cluster in a user-centric way [21-23]. Each 
user reports a suggested CoMP cell set based on the reference 
signal receive power (RSRP) measurement from each cell, where 
a non-serving cell is included in the reported CoMP set only if 
its corresponding RSRP is within the given offset of the RSRP of 
the serving cell, represented by [21] 

Cell i G reporting CoMP set, if and only if 
RSRP, > RSRPo - a th , 

where RSRP denotes the RSRP of the serving cell, and RSRP £ 
denotes the RSRP of the cooperating cell reported by the user, 
and ath denotes a given threshold whose value is usually set to 
a small portion of RSRPo. Then, CoMP user scheduling can be 
performed and CoMP set will be determined correspondingly. In 
reality, it is possible that different users report overlapping CoMP 
sets [21-23], in this case iterative scheduling might be used to 
achieve a desirable performance. If ath is small, it is understood 
that the scheduled CoMP users see significant interference from 
the cooperating cells, where a coordinating transmission is highly 
desirable. Here we focus on the design of coordinated multicell 
beamforming within a single cell cluster. Consider a coordinated 
cell cluster composed of K users and B BSs. Each BS is equipped 
with M transmit antennas and each user is equipped with a single 
antenna. Assume that each BS has at least one served user, i.e., 
K > B. The transmit signal from each BS is only intended for 
its own users, meaning that no user-data sharing between the BSs 
is needed. However, the transmit beamformers from different BSs 
need to be designed jointly to suppress the ICI. 

Let ICb denote the user set served by BS b, and denote the 
information symbol intended for user k, with E{|sfc| 2 } = 1, where 
E{ } denotes the expectation operator. Let bk be the index of the 
BS serving user k, ffc G C Mxl be the unit-norm beamforming 
vector (i.e., ||ffcj| = 1 Vfc) for user k and pk be the transmit power 
for user k. Each BS is restrained by a power constraint Pi,, i.e., 
J2keK P k — ^ BS - Assuming k e ICb, the received signal of user 
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k is given by 
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where f denotes Hermitian transpose, n k is the zero-mean com- 
plex white Gaussian noise with variance a\, and hfc& denotes 
the frequency-flat fading channel vector from BS b to user k. 

,Pk] T with the superscript T denoting the 
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and W = diag(fi,f 2 , . . . , f^)- The 
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where U(j) denotes the index of the BS serving user j, Note that 
although the above model appears to assume perfect symbol-to- 
symbol synchronization, this assumption can be removed in or- 
thogonal frequency-division multiplexing (OFDM) systems when 
the SINR or some utility function monotonically increasing with 
SINR is used as the performance metric because the SINR 
criterion only involves the power of the interference. 

The problem considered in this paper, is the joint optimization 
of the BS beamforming vectors {f k }y k and the power allocation 
{Pk}\/k- A TDD multicell system is first investigated, where 
each BS (say b) can directly estimate the CSI of its own users, 
i.e., {h. k b} k <=K. b , without information exchange between cells, by 
exploiting downlink-uplink reciprocity. It is also assumed that 
each BS can estimate from the reverse link the crosstalk channels 
to other users in the same cluster. This is justifiable because 
usually the threshold ath is small and therefore the channel from 
the BS to the users in other cells and that from the BS to its own 
users are of similar order of magnitude. To do so, a coordinated 
training phase for the cell cluster might be required, in which the 
training signals from different users are designed to be orthogonal 
in order to avoid mutual interference. In summary, {^i k b} k= ,i is 
referred to as local CSI at a TDD BS b and will be exploited 
in the subsequent design. Note that our focus is not limited to 
TDD multi-cell systems. Our developing schemes can be easily 
generalized to frequency division-duplex (FDD) systems via CSI 
feedback from user terminals, to be discussed later. 

III. Centralized Multicell Beamforming 

Various criteria have been used to design multi-cell beamform- 
ing schemes, including the sum-rate maximization, the worst-user- 
rate maximization and the sum-power minimization with SINR 
constraints. In contrast to these objectives, our design focuses 
on the tradeoff between rate efficiency and fairness, as each 
coordinating cell would expect that a certain cooperation gain is 
guaranteed. If we employ a game-theoretic model and view users 
in the multi-cell coordinating systems as players with conflicting 



objectives, it can be shown that the cooperative NBS could achieve 
a reasonable compromise between efficiency and fairness [10, 24, 
25]. However, the algorithms usually require sharing of global 
information between the BSs and it is hard to realize a fully 
distributed solution. With the goal to achieve distributed multicell 
beamforming optimization, we alternatively resort to the max-min 
fairness based optimization and take the objective as the max-min 
point on the Pareto boundary of user rates. To gain some insight 
on our distributed solution, we will first provide a centralized 
algorithm consisting of two steps. The first step solves a weighted 
SINR balancing problem in order to reach the max-min fairness 
with the least sum-power, which is followed by the second step 
of Pareto improvement targeting at the Pareto optimality. The 
detailed algorithm is described as follows. 

A. Two-step Optimization 

To achieve the max-min fairness, here we take the minimal 
SINR as the objective function, i.e., to maximize the worst-user 
SINR. The optimization problem is formulated as 



^SINR^p) st Epfe < PBsV6i (5) 

Pk k£K b 



W,p 



where p k is a weighting factor to take into account heterogeneous 
requirements from users; users with larger factors will achieve 
proportionally higher SINRs. Note that though here we adopt the 
minimum weighted SINR as the objective function of the max-min 
optimization, i.e., to maximize the worst-user weighted SINR, this 
is equivalent to the max-min rate optimization, as the transmission 
rate is monotonically increasing with the SINR. Problem (5) is 
difficult to solve due to its non-convexity. Therefore, we first ex- 
plore its relationship with the following sum-power minimization 
problem 



SINR? L 



mm||p||i s.t. 
w ,p 



> 7 Vfc, 
Pk 



(6) 



in which 7 is a preset target and || • ||i returns the 1-norm of 
an input vector. It is easily seen that the solution to (6) attains a 
set of equal weighted SINRs. If we denote the optimal weighted 
SINR of (5) as 7* and the corresponding beamforming vectors and 
power vector as W* and p*, respectively, then it can be proved, 
see [26], that the solution to (6) with a target 7* will produce 
the same optimal W* and p*. As a result, if (6) is solved, (5) 
can be simultaneously solved by searching for the highest target 
7* for which the solution p* of (6) has at least one active power 
constraint. Furthermore, as shown in [20], (6) can be transformed 
into a standard second-order cone programming (SOCP) problem 
and hence can be solved using standard optimization packages, 
such as CVX [27]. Based on this, a bisection search method can 
be employed to efficiently search for the optimal weighted SINR 
target 7*, hence giving the solution to (5). 

Numerical results show that to achieve the above solution, in 
many cases only a subset of coordinating BSs need to transmit at 
full power, while other BSs transmit with less power. Therefore in 
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general it is possible to further improve the sum-rate without los- 
ing the max-min optimality. This implies that the above weighted 
SINR balancing solution is not Pareto optimal. Motivated by this, 
we next propose a Pareto improvement method in order to reach 
the Pareto boundary with max-min fairness. 

Assuming M > K, the proposed Pareto improvement is 
realized by the following beamforming updating algorithm. With 
a slight notation abuse, let {f/!}vfc and p* = \p\, . . . ,p* K \ T 
denote the unit-norm beamforming and power solutions to (5), 
respectively. For BS b, if J2ieK b Pi < P bs and p* k = min ie!Cb p*, 
then we propose that user k transmits at full power and its unit- 
norm beamforming vector is updated as 



r new 
I k 



(7) 



method achieves the maximum rate improvement, but there might 
still be a possibility that an update along some other subspace 
gives a further rate improvement. Though the proposed algorithm 
cannot ensure Pareto optimality in the general case, the following 
theorem shows that Pareto optimality with max-min fairness 
indeed is achieved in the two-user case. 

Theorem 1: For the two-BS two-user coordinating system with 
M > K = B = 2, the proposed two-step optimization algorithm 
guarantees the Pareto optimality with max-min fairness. 

Proof: See Appendix A. ■ 
Note that if the proposed weighted SINR balancing optimization 
suggests that one BS transmits at less power, the proposed 
beamforming update method would achieve a maximum rate 
improvement and thus guarantee the Pareto optimality. 



where 8 k = Z(hLf£), a k is a positive scalar such that ||f / new || = 
1, and h| F is the projection of h k t, onto the complement of the 
column space of H fe = [hi 6 • ■ ■ h k _ lyb h fe+ i, 6 • ■ ■ h Kb ], given by 



H fc (H[H fc 



si 



(8) 



This completes the proposed two-step multicell beamforming 
scheme. 

B. Pareto Optimality 

Here we will discuss the effectiveness of the above proposed 
scheme, both in terms of max-min fairness and Pareto optimal- 
ity. In the first step the weighted SINR balancing optimization 
guarantees the max-min fairness, weighted by the factors {p k }. 
We will prove that in the second step the Pareto improvement 
method yields a rate increase without losing the achieved max- 
min fairness. For the special case of two BSs and two users, 
i.e., B = K = 2, Pareto optimality of the proposed scheme is 
guaranteed. 

To proceed, we first look at the interference of the signal of user 
k on other users. The property h| fc h| F = yields that, with the 
updated beamforming vector f£ ew , the interference power caused 
on any other user i ^ k is given by 

2 
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P k 



Vi T f* 
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This means that the beamforming update in the second step will 
not change the interference level imposed on other users, i.e., it 
will not harm the max-min optimality achieved by the weighted 
SINR balancing. However, the beamforming update is able to 
further improve the user's rate. The effective power of user k with 
the beamforming update can be written as (10), shown at the top 
of the next page, where equality holds only if the two terms on 
the right-hand-side have the same phase, i.e., 9 k = /(h^f^). 
That is, the proposed beamforming update method is able to 
improve the rate as much as possible while keeping the max- 
min optimality. To further illustrate the efficiency of the proposed 
two-step optimization algorithm, it is interesting to investigate its 
Pareto optimality. 

As seen from (7), the proposed beamforming update for user 
k £ JCb restricts its beamforming update only in the null subspace 
of Hfc. Under this subspace restriction, the proposed update 



IV. Distributed Multicell Beamforming 

The above solution is by itself centralized as far as implemen- 
tation is concerned. Here, our aim is to develop a distributed 
multicell optimization algorithm that reduces the amount of in- 
tercell communication. For the step of Pareto improvement, it 
is seen from (7) that the beamforming vector update for each 
BS can be performed distributively using local CSI. The main 
difficulty, however, lies in the distributed implementation of the 
weighted SINR balancing optimization, which will be addressed 
by investigating an approximate uplink-downlink duality. 



A. Uplink-Downlink SINR Duality 

If we view the cluster of BSs as a super BS (SBS), our system 
model can be regarded as a virtual multiuser system consisting of 
K user terminals and a SBS. 2 Previous results in [28,29] have 
shown that for each user, say user k, the same SINR target can 
be achieved in both downlink and uplink with the same set of 
beamforming vectors and the same sum-power. This implies that 
only if a sum-power constraint for the virtual multiuser system is 
considered, the existing uplink-downlink duality can be directly 
applied to solve the downlink optimization problem by converting 
into a virtual uplink problem, which enables a distributed imple- 
mentation. When per-BS power constraints are considered, the 
conventional uplink-downlink duality becomes inapplicable. It is 
worth mentioning that a modified form of duality proposed in 
[28] is able to deal with the issue of per-BS power constraints 
and find the exact optimal solution. However, the subgradient- 
based algorithm in [28] shows a very slow convergence and this 
brings difficulties to implementation. As shown in [17], the inter- 
BS communication required in its implementation involves about 
O(NK) positive scalars, where the required number of iterations 
TV is typically about one hundred as shown, thus imposing a heavy 
burden on the backhaul link. Balancing the pros and cons, here, we 
propose to use an approximate uplink-downlink duality to achieve 
a distributed multicell beamforming scheme with limited intercell 
communication. 



It should be noted that data sharing between the BSs is not permitted in our 
system model. Therefore, the multiuser precoder of the SBS should be restricted 
to be in a block-diagonal form expressed as diag{fi, . . . , fx}- 
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With the goal to solve the max-min optimization, we develop a 
virtual uplink problem which is approximately dual to (6) and fa- 
cilitates a limited-intercell-communication distributed realization, 
given by 



mm q i s.t. 



SINR ^ L 

pk 



> 7 Vfc, 



X) 8fc < Pes Vh, 



(11) 



fce/c b 



in which q = [qi, . . . , qK\ T denotes the virtual uplink power 
vector and 



SINR^ L (f fc ,q) = 



E <liH-i 



I ft 



(12) 



denotes the uplink SINR of user fc. Next we will clarify the 
approximate duality between (11) and (6). 

Denote the solution to (11) as W* and q* = [q*, . . . , q* K ] T , 
where q* satisfies the per-cell power constraints Efce/Ct, 9fc — 
Pbs- We convert this uplink power solution q* to the downlink 
using the uplink-to-downlink converting method in [29], and 
denote the downlink power solution as p* = \p*, , . , ,p* K ] T , which 
satisfies ||p*||i = ||q*||i- The conventional uplink-downlink dual- 
ity theory reveals that the uplink SINR achieved by W* and q* 
is equal to the downlink SINR achieved by W* and p*. Note that 
strictly speaking W* and p* are not exactly the solution to (6), as 
p* cannot be guaranteed to satisfy the per-BS power constraints. 
However, it is found that for any achievable 7, the equality 

max & J2keK b Pl = max bEfceKi, Ik ^ p BS approximately holds 
in most cases, as shown in Fig. 3 (which will be mathematically 
justified later), implying that p* approximately satisfies the per- 
BS power constraints YkeK b Pk — Pbs- This property reveals 
an approximate duality between the uplink problem (11) and the 
downlink problem (6), i.e., the same SINR can be achieved by 
the uplink and the downlink with the same set of beamforming 
vectors and approximately the same individual (per-cell or per-BS) 
power constraints. Based on this observation, the solution to the 
downlink problem (6) can be approximately achieved by solving 
the uplink counterpart (11). 

To mathematically justify this approximate duality, we define 



D(W*) = diag • 



PKl 



(f*) 7ii,n(i)f* (f^ H k ,u(k)^k 



(13) 

and a K x K matrix Vl^W*) with its (i, j)-th element given by 



(f;)t^ jU0) f; for i f j, 

for i = j. 



(14) 



It was shown in [29] that the downlink and uplink power vectors 
can be obtained as 



P 

q* 



Al, 
A T 1, 



(15) 
(16) 



where 



(17) 



A = [I - D(W*)*(W*)] _1 D(W*) 

= [D _1 (W*) - *(W*)] _1 , 

and 1 is the all-one vector. It is difficult to derive a rigorous 
condition for the uplink-downlink duality with individual power 
constraints, i.e., max 6 YkelC b P*k = max b EkeK b € < p bs- As 
an alternative, our analysis resorts to examining how far from 
being diagonal A = D _1 (W*) - *(W*) is, as the duality holds 
exactly when it is strictly diagonal. That is, we will examine the 
duality between (6) and (11) by seeing if A is diagonal. We will 
show that (11) is diagonal-dominant or quasi-diagonal in most 
cases, i.e., the duality approximately holds. 
To proceed, we define the following metric: 

(f *\> U -F* 
(IfeJ «fc,U(fe)I/c 



Vk 



(18) 



which is the ratio of the magnitude of diagonal entry in the k- 
th column to the sum of the magnitudes of other entries in this 
column. A matrix A is said to be strictly diagonal-dominant if 
r] k > 1 Vfc [30]. Recalling that {f£} is the solution to (11), we 
therefore have 



Pkl 



< 



Based on the above, rjk can be rewritten as 



I f,* 



> 



i,n(fe)ffe 



(19) 



(20) 



Note that typical user-centric cell clustering strategies [21-23] 
result in that user-BS links have comparable channel strengths, 
resulting in that {q*}f =l will have the same order of magnitudes 
to yield rj k > 1 Vfc. 

It is also useful to examine the asymptotic behavior of A. When 
the received SNR becomes extremely small, i.e., ^'Hi.u(fc) = 
<? fc h, iI]( fc)h. iI](fc) ^ ^ ^ ^ ^ ^ ?y fe > 1. This 

implies that A and also A are quasi-diagonal. Similarly, when the 
received SNR becomes extremely large, the inter-user interference 
becomes dominant over the noise. Thus, the main purpose of the 
beamforming optimization in solving (11) now lies in mitigating 



the inter-user interference, leading to (f^) Hi^(k)^k 
This yields again rj k >• 1. 



OWi^k. 
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B. Distributed Implementation 

To achieve the weighted SINR balancing problem (5) in a 
distributed way, here we propose an iterative algorithm to achieve 
the solution to the approximate virtual uplink problem, and then 
convert it to the downlink. Based on the derived approximate du- 
ality, the virtual uplink problem corresponding to (5) is formulated 
as 



. SINR^ L (W,q) 

max mm s. 

w, q k p k 



BS 



V6. 



(21) 



Using a similar method as the one in solving (6) in Section 
III- A, (21) can be solved by iteratively solving the sum-power 
minimization (11) and searching for a maximum achievable SINR 
target. As shown in [17], such a strategy permits a distributed 
solution using a standard power control loop [31]. However, its 
implementation requires a two-layer iteration and suffers from 
slow convergence. To speed up the convergence, we propose a 
single-layer iterative algorithm to reach a distributed solution. 
Proposed Iterative Algorithm 

0) Use the superscript (n) to denote the corresponding param- 
eters in the rtth iteration. 

1) Initialize q^ - 1 with drawn from a uniformly distributed 
random variable between (0, Pbs] or a simple all-zero vector, 

v (0) = 10 -3 and n = l_ 



T 

2) Set S < n) 
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"^//.(q^™ for all k, where the interference 



function Ik{q^ n ') is defined as 



4(q (n) ) 



in which 



life ) ] rik,n(k)h 
9fe n)f fc^fe,u(fc) f fc 



are; max — , 

'•• n (Eh <i 



n 



i,n(fc) 



I ffc 



(22) 



(23) 



3) Find a = mirif, 



(n) 



4) Update q y k 

5) Find the 



= ag^ for all k. 
minimum SINR 



and 



set 



min fc SINR^f^W"' 



6) Stop if the convergence condition defined by Hq*-™-* 



< e (a positive scalar close to zero) is satisfied. 



Otherwise, n = n + 1 and go back to Step 2. 

Theorem 2: The proposed iterative algorithm with any arbitrary 
initialization always converges to the global optimum of (21). 

Proof: We define 1(7, q) = [7/1 (q), . . . ,7/if(q)] T , which 
can be recognized as a standard interference function [31] for 
the SINR target 7. We also find it useful to first review the 
key properties of the fixed-point iteration algorithm. Assuming 
a feasible initial point q(°' to achieve 7, the fixed-point iteration 
q(n) _ x(7,q(" -1 )) has the following properties [31]: 

PI: q(") is component-wise monotonically decreasing; 
P2: q(") converges to the unique optimal solution; 
P3: qW for all n are feasible solutions. 

Our proof is divided into two parts. We first prove its con- 
vergence to a fixed point, and then prove that this fixed point 



is optimal. We start with Step 5 in the nth iteration. Obviously, 
q(") satisfies the individual power constraints and achieves 7W. 
At Step 2 in the (n + l)th iteration, g("+ x ) = X( 7 ("), q(™)) is 
obtained via a standard mapping. Then, according to PI, we have 

(24) 



r("+l) < q («)_ 



As a result, at the following Step 3, a > 1. Based on the SINR 
expression in (12), this yields 

.UL^(«+l) „, .(n+X)' l 

(25) 



SINR^ L (f^ +1) ,q("+ 1 )) = SINR£ L (f^,ag^ +i 
>SINR^ L (ff +1) ,g(" +1 )) 



According to Step 2 and the definition of we have 

mmSINRU L (f(" +1) ,g("+ 1 )) > 7 H (26) 

k 

Therefore, ^ n+ ^ > i.e., the minimum SINR is monoton- 
ically increasing and the algorithm converges to a fixed point 
{7^°°', q(°°'}, in which at least one individual power constraint 
should be active. Next, we proceed to prove that the fixed point 
is also optimal by contradiction. 

Assume that the optimal solution is {yP^q ?*} and -yC 00 ) < 
7° pt . We see that q^°°^ satisfies that at least one individual power 
constraint should be active and 

g («0 = 7 Mj fc ( q (°°)) Vfc. (27) 

As a consequence, q(°°) is the unique solution to (11) with the 
SINR constraint according to P2. Now we take q opt , which 
is also a feasible solution as the initial vector for the fixed-point 
algorithm to solve (11) and the converged power vector is q. 
Then, according to PI, we have q < q opt which is expected 
to be the same as q^°°\ Since q opt should satisfy the individual 
power constraints, this indicates that with q, no individual power 
constraint is active and therefore contradicts that q(°°) is the 
unique optimal solution to (1 1). ■ 
After we obtain the solution of (21), say {7, q, W}, the solution 
to the downlink problem can be achieved using a standard power 
control loop in a distributed manner. However, in order to fully 
utilize the backhaul link between the cooperating BSs, here we 
choose to calculate the downlink power vector based on the 
following linear system of equations [20, (19)— (21)] 

p = [D _1 (W) - *(W)] _1 1. (28) 

It is possible that some pk obtained from above could violate the 
individual BS power constraints, due to the fact that the duality 
with individual power constraints only holds approximately but 
not exactly. In that case, the power of the specific BS is scaled 
down to fulfill the power constraints, which would possibly result 
in some performance loss in terms of the max-min fairness. 
Our simulation results in Section V will demonstrate that the 
percentage of the power constrain violation is marginal and its 
resulting worst-user rate performance loss is insignificant. On 
the other hand, it is worth mentioning that the possible scaling 
down of the power solution to the max-min optimization problem 
would provide more room for the sum-rate improvement in the 
second step. Simulation results show that our distributed scheme 
usually achieves a better sum-rate performance though suffering 
a performance loss in terms of the worst-user rate, compared to 
the centralized algorithm. 



7 



C. Required Inter-BS Communication 

In this subsection, we analyze the required inter-BS communi- 
cation in the proposed distributed algorithm. As has been noted 
before, the Pareto improvement can be separately implemented at 
each BS using local CSI {hfc&}£L-,, without the need of inter-BS 
communication. The information exchange between the BSs is 
only required in the max-min optimization in which the proposed 
iterative algorithm is first implemented and its solution is then 
converted to the downlink. 

Assume that there are B coordinating BSs simultaneously 
serving K users, each of which (say BS b) has \ICb\ users. By 
examining the steps of the proposed iterative algorithm, we can see 
that the results will not be changed if we use — /^(q'"^ 1 )) in 
Step 2 instead of g { k n) = h^ 1 - 1 ')), so is not needed 

in actual calculation. As a consequence, the implementation only 
requires that the power vector q( n ) at each iteration be shared 
between the BSs. That is, each iteration requires the exchange of 
Ylb=i (B— l)l^-f>l = K(B— 1) positive scalars between the BSs 3 . 
Suppose that the algorithm converges within (N + 1) iterations 4 , 
the inter-BS communication involves NK(B-l) positive scalars. 
Then we look at the computation of the downlink power vector p 
from the converged result of the iterative algorithm, it can be seen 
from (28) that each BS can distributively calculate its preferred 
power by sharing matrices D and \&. As shown in the definitions 
of D and it needs K(B - 1) ^f =1 \JC b \ =K 2 (B- 1) positive 
scalars of inter-BS communication. 

To highlight the light inter-BS communication requirement of 
the proposed method, we compare it with a simple strategy where 
full CSI is shared between the BSs and each BS individually 
performs beamforming optimization using global CSI information. 
Recognizing that exchange of one M-dimension channel vector 
requires about AM positive scalars, the latter strategy would take 
about AMKB(B — l) positive scalars of inter-BS communication, 
which increases linearly with the number of antennas at the BSs. 
In comparison, the ratio of the inter-BS communication of the pro- 
posed method to that of full CSI exchange is calculated as tmW - 
Our simulations show that the convergence speed of the proposed 
method is not significantly affected by the number of BS antennas, 
due to the fact that the number of the dual variables in our 
problem formulation only depends on the number of coordinating 
BSs and users but not the number of antennas. This reveals that 
our proposed method has much reduced requirement for inter-BS 
communication, especially when the number of transmit antennas 
is large. 

D. Practical Issues 

Note that the proposed distributed scheme assumes that each 
BS b knows its local CSI {hjy,} vi . perfectly. This is hard to 
realize in practical systems. It is more reasonable to assume that 
only imperfect CSI is available at the BSs. In particular, if the 
bandwidth of the backhaul link is limited, the implementation of 

3 This is evaluated based on the signaling overhead of backhaul links. If one 
BS needs to transmit one positive scalar to other 73 — 1 BSs, it would involve an 
exchange of B — 1 positive scalars. 

4 Our simulation results show that typically around 5 iterations are sufficient to 
achieve a desired accuracy. 



the proposed iterative optimization will suffer a fairly large delay 
because of the parameter exchange between the BSs. Thus, the 
CSI could be outdated by the time the iterations are completed, 
resulting in some performance loss. This problem becomes severe 
in high-mobility scenarios where the channels vary rapidly. Even 
in low-mobility scenarios, it is still possible that CSI varies while 
the algorithm is running. Moreover, a certain level of channel 
estimation error is inevitable in practice. This motivates us to 
evaluate the robustness of the proposed scheme against imperfect 
CSI. 

Imperfect CSI is commonly modeled either stochastically or 
deterministically. In the former case one assumes that the un- 
certainty region of the CSI errors is unbounded and distributed 
according to some known distribution [32-34], while in the latter 
case the uncertainty region of the CSI perturbations is assumed to 
be bounded [35,36]. Here we will adopt the bounded imperfect 
CSI model, as we can evaluate the performance of our proposed 
scheme under a restricted iteration delay. This can be realized by 
fixing the number of iterations, such that the backhaul signaling 
overhead is restrained and the resulting delay is also limited. 
This kind of evaluation is of particular importance for practical 
applications. By decomposing the channel vector into a large-scale 
factor and a small-scale fading part, the erroneous channel can be 
expressed as 

h kb = Ttohjg = lkb (hjg + h«) , (29) 

where jf-j, denotes the large-scale channel component and is 
assumed to be exactly estimated by the BS, hj^ denotes the actual 
small-scale fading part of the channel, denotes the nominal 
small-scale fading channel known by the BS, and h^ 1 denotes 
the channel uncertainty. As is common in the literature [35], we 
assume that the channel uncertainty is confined within an origin- 
centered hyper-spherical region of radius a^b, i.e., ||h!X||2 < &k b - 
Simulation results in Fig. 8 will show that for a CSI uncertainty 
with o~kb = 0.1, the proposed distributed solution for M = 4 
only suffers from marginal performance loss. It is also worth 
noting that, the strategy of robust design in [35] can be readily 
extended to our scheme. To do so, we only need to define 
sinri! 1 = minn.rfs),, , , SINRF L as the worst-case SINR over 

the uncertainty regions, and use it to substitute SINR^ L in the 
objective function. Several possible bounds of sinr)! 1 , e.g., [35, 
(13)], can be used in our algorithm. 

Finally, we will discuss the implementation of our scheme in 
FDD multicell systems. Different from TDD implementation, in 
FDD mode the CSI can only be obtained through feedback from 
users. It is noted that normally each user terminal, say k, can only 
estimate its channels from all the coordinating BSs, i.e., {hk b } vb , 
and feed back these CSI to its serving BS. Therefore, in order to 
generalize our scheme to FDD systems, extra channel information 
exchanging between the BS will be needed such that each BS b 
can have CSI {h-kb}yk and then perform iterative beamforming 
optimization. 

V. Simulations Results 

In this section, computer simulations are provided to evalu- 
ate the performance of the proposed schemes. The user-centric 
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cell clustering strategy described in Section II is used in our 
simulations. In particular, our user scheduling only considers the 
cell-edge users reporting the same 2-BS or 3-BS CoMP set with 
a>th = 0.4 RSRPq. More specifically, we assume that the inter-BS 
distance is 1km and the users have at least 350m distance from 
their serving BS. For each time slot, each BS uses a round-robin 
strategy to schedule a single user from its serving user pool. We 
consider fiat fading channels, where the small-scale fading is 
assumed to be zero-mean Gaussian distributed with covariance I, 



i.e., h 



(s) 
kb 



CAf (0,1), the large-scale path loss jkb is given by 



0/i 



hi, 



(30) 



in which j3 = 10 3 45 is a scaling factor, and I = 3.8 denotes the 
path loss exponent, dkb is the distance between user k and BS b. 
In dB, this gives 10 log 10 (7 fc6 ) = -381og 10 (d fc6 ) - 34.5 + (i kb , 
where the shadow fading follows a normal distribution ^b ~ 
Af(0, 8 dB). Throughout the simulations, we set p k = 1 and 
assume that the coordinating BSs have the same power constraint, 
and the scheduled users have the same noise figure. The user SNR 
is defined as the ratio of the average signal power from its serving 
BS to the noise power. For simplicity, we assume the scheduled 
users have the same path-loss to their serving BSs and thus have 
the same SNR. 

For comparison, several existing multicell beamforming 
schemes based on different criteria are also simulated, including 
the SGINR [18, 19] scheme, the sum-rate optimal scheme and 
the NBS. The optimal sum-rate is realized by numerical search 
over the achievable rate set generated using the method in [11, 
Corollary 2]. The NBS is achieved by solving the optimization 
problem as follows [10,24,25] 



k 



s.t. 



J2 Pk < Pbs Vb, (31) 



where {R®} denotes the initial agreement point. As is commonly 
done, we use the Nash Equilibrium (NE) solution {-R^ } as the 
initial point. The above NBS problem in general is non-convex, so 
we employ alternating optimization between W and p to find a 
near-optimal solution. Taking into account the possible limitation 
of backhaul capacity in practical systems, we also simulate the 
proposed scheme with a fixed number of iterations, such that the 
information exchange between the BSs is limited. 

The Pareto boundaries of the multicell beamforming system 
with (M, K) = (4, 2) under two random channel realizations are 
first simulated and shown in Figs. 1 and 2, using the approach [11, 
Corollary 2]. It is seen that except the NE point, other schemes 
all achieve Pareto optimality in these two cases. In particular, 
as shown in Fig. 1, the rate point achieved by the proposed 
centralized scheme is located exactly at the intersection point 
between the Pareto boundary and the line R2 = R\. In another 
case as shown in Fig. 2, the Pareto boundary has no intersection 
point with the line R2 = R\, implying that the SINR balanced 
rate point does not operate on the Pareto boundary, but the rate 
point of our centralized scheme is still on the Pareto boundary 
and is the point closest to the line R 2 = R\. This means that in 



these two cases the rate point achieved by the proposed centralized 
scheme always has its worst rate not less than any other points 
on the Pareto boundary, guaranteeing a strict max-min fairness. 
As for the proposed distributed scheme, it is seen that the rate 
point also is located on the Pareto boundary, between the rate 
point of the centralized and the sum-rate optimal point, having a 
behavior similar to the NBS solution. This implies that compared 
to the centralized scheme, the distributed scheme has a certain 
performance loss in terms of max-min fairness, due to the fact 
that the duality is only approximate, although it achieves a better 
sum-rate performance. 

The analysis in Section IV-A has shown that if the condition 
maxf, X)feeAC t Pk — -^bs is met , the optimization of the virtual up- 
link problem can guarantee the efficiency of the resulting downlink 
solution. Otherwise, if max ^ exceeds the individual power 

b k£K b 

constraint Pes, our distributed scheme will in general suffer from 
a performance loss in terms of max-min fairness, depending on 
how much the power constraints are exceeded. For this reason, 
we define the following excessive power percentage 



max b Efc GKb P*k - A 



BS 



Pbs 



(32) 



Fig. 3 shows the cumulative density function (CDF) of tp under 
several configurations. We see that max^-p^. < Pbs holds with a 
probability of more than 70%, and <p has a value below 3% with 
a probability of more than 90%, showing that the duality holds 
with high probability. 

Fig. 4 illustrates the average rate results achieved by the worst 
user in the proposed multicell beamforming schemes and the 
reference schemes for (M, K) — (4, 3), calculated over 5000 
random channel realizations. The results reveal that the proposed 
centralized and distributed schemes both outperform the reference 
schemes over a wide SNR range. In particular, the centralized 
scheme achieves a gain around 0.5 bpsk/Hz over the sum-rate 
optimal scheme, the SGINR scheme and the NBS. In contrast 
to the centralized scheme, the distributed scheme shows a per- 
formance loss around 0.25 bpsk/Hz. Results also show that most 
performance improvements in our distributed scheme are achieved 
already within 2 iterations. This means that the efficiency of our 
distributed scheme can be maintained even with very limited 
inter-BS communication. Fig. 5 shows the CDF of the worst- 
user rate results achieved by the above schemes. It is seen that at 
all outage levels, our centralized scheme exhibits the best worst- 
user rate performance. Though our distributed scheme suffers a 
slight worst-user rate loss compared to the centralized scheme, it 
outperforms the NBS and the sum-rate optimal scheme. Looking 
at the sum-rate performance, Fig. 6 provides the CDF of the ratio 
of the sum-rate achieved by the above schemes to the optimal 
solution. The results demonstrate that the guarantee of strict max- 
min fairness by our centralized scheme comes at the cost of 
a large sum-rate performance loss. On the contrary, our two- 
step optimized distributed scheme achieves a significant sum-rate 
improvement over the centralized scheme by sacrificing a certain 
max-min fairness. It is seen that the outage behavior of the sum- 
rate achieved by our distributed scheme is very close to the NBS. 
At an outage probability of 10%, more than 80% of the optimal 
sum-rate is achieved. Combining the results in Fig. 5, in this setup 
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a random channel realization, with (M, K) = (4, 2) and SNR=15 dB 



a random channel realization, with (M, K) = (4, 2) and SNR=0 dB. 



our distributed scheme outperforms the NBS in terms of fairness 
rate and has a comparable sum-rate performance. 

To evaluate the impact of quantized parameter exchange be- 
tween the BSs on the performance of the proposed distributed 
scheme, we provide the average worst-user rate performance of the 
distributed scheme applying a simple linear uniform quantization 
between and Pes to the power parameters c£ n \ It is seen in 
Fig. 7 that a 3-bit quantization gives an inappreciable performance 
loss. Finally, to investigate the robustness of our distributed 
scheme against the CSI imperfection, we adopt the CSI error 
model (29) and set M = A,a kb = {0.05,0.1,0.2}, plotting 
the average rate of the worst case user assuming that the CSI 
errors are uniformly distributed in the uncertain region. S imulation 
results in Fig. 8 show that though CSI errors generally degrade the 
performance and the loss increases with the SNR, no error floor is 
observed at SNR = 20 dB in any of the setups. A CSI error level 
less than <Jkb = 0.1 gives a slight performance loss especially 
in low and medium SNR regimes. Even with <Jkb = 0.2, the 
proposed distributed scheme still outperforms the SGINR scheme 
with accurate CSI over a wide SNR range. 

VI. Conclusion 

This paper has presented two fair-rate optimal multicell down- 
link coordinated beamforming schemes. We first devised a two- 
step centralized optimization scheme, in which the first step 
achieves a strict max-min rate fairness, while the second step 
further improves the sum-rate as much as possible without losing 
the max-min optimality. We proved that this scheme is guaranteed 
to reach the Pareto boundary with max-min fairness in the case of 
two BSs each serving a single user. By exploiting the approximate 
uplink-downlink SINR duality with individual power constraints, 
a distributed solution, in which only a small number of positive 
scalars are shared between cells, was proposed. Simulation results 
showed that our distributed scheme can approach the max-min 
Pareto optimality, requiring only very limited inter-BS communi- 
cation, suitable for practical BS cooperation. 
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Fig. 7. The average worst-user rate of the proposed distributed multicell scheme 
employing uniform linear quantization on the exchanged power parameters, and 
fixing 2 iterations, K = 3. 



Appendix A 
Proof of Theorem 1 

We consider two separate cases. If the solution to the weighted 
SINR balancing problem suggests that both BSs use full transmit 
power, then it is impossible to obtain a Pareto improvement, 
meaning that the Pareto optimality with max-min fairness has 
already been achieved. Therefore, we only need to focus on the 
other case where one BS uses full transmit power while the other 
one uses less power. Without loss of generality, assume p\ < Pbs 
and p* 2 = Pbs ■ The proof of Theorem 1 starts with the following 
Lemma. 

Lemma 1: In the solution to the weighted SINR balancing 
problem with M > K = B = 2, suppose < Pes and 
P2 = Pbs- If a Pareto improvement over the balanced rate tuple 
(r*,r*) can be achieved in which user l's rate improves while 
user 2's rate remains, then the updated beamforming vector for 
user 1, denoted by ff cw , must generate no interference to user 2, 
i.e., f 1 nowt ?i 21 f 1 now = 0. 

Proof: The proof is by contradiction. First, we let the 
beamforming vectors and the transmit power of user 2 be fixed by 
f" cw , f| and Pbs, respectively. Then, we can express the SINRs 
(assuming pk = 1 for simplicity) of the two users as functions of 
Pi, given by 

Pill "-lltl 



SINR? L ( Pl ) 



SINR° L ( Pl ) 



Pbs^H-i^ 
p*t 



1 



^BSh fi -22h 



ftfr'^ifr + 1 



(33) 



(34) 



Define ASINR(pi) = SINR° L ( Pl ) - SINR^"" (pi). We observe 
that ASINR (p\ ) is a monotonically increasing and continuous 
function. The assumption of Pareto improvement thus yields 



>DL 



ASINR (p" cw ) > 



but it can be easily observed that ASINR(O) 
consequence, there must exist some p[ < 
< ASINR (pi) < ASINR (p" GW ). Now, if we assume that 



(35) 

< 0. As a 
pi lcw such that 



II 



f«ewt^ 2if new ^ q is true, then this implies that if BS 1 chooses 
to use transmit power p[, the resulting SINRs of both users will 
be greater than 7*, which contradicts the result that 7* is the 
solution to the max-min optimization. Therefore, we must have 
f^ ew '"H 2 iff cw = 0, which completes the proof of Lemma 1. ■ 

Now we proceed to prove Theorem 1. Lemma 1 reveals 
that the updated beamforming vector for user 1 should satisfy 
£newt^£ 2i fncw _ q -p n j s a j so implies that it is impossible for the 
SINR of user 2 to improve over 7* and any interference to user 
1 will lead to the result that its achieved SINR becomes smaller 
than 7*. As a by-product, it can be seen that f* T^if* = 0. 

Based on these observations, if we denote the beamforming 
update as 

Afj - VPp f r w - y/pft, (36) 

we must have h^Afi = 0. Further, it was shown in [11] that 
the Pareto optimal beamforming vector for any user always can 
be expressed as a linear combination of the channels between its 
serving BS and all the users. Following a similar procedure, it can 
be proved that the max-min optimal beamforming vector f£ has 
the same property. That is, they both can be formulated as 



fl = C2lh; 

C21 + Cl 



,ZF 



21 -I- 1,111111 = (,21«21 -I- t,ll«l , (37) 

where C21 = C21 + Cn^fef • It follows from fj%ifi = that 
C21 = 0. This means that the beamforming update should be 
restricted in the direction of hf F . Thus, the updated unit-norm 
beamforming vector f" ow and the updated transmit power p" ew 
can be written as 



yPi r i 



(38) 



where b\ is a complex scalar. For the maximum Pareto improve- 
ment, the optimal power allocation p ncw and the scaling factor 61 
can be determined by solving 

2 

f new 
r l 



0<pj° 



max 

b t ec 



VPl 



ncwu I c 11 < 
n ll r l 



S.t. 



1. 



(39) 



Based on the triangle inequality, it can be readily proved that the 
proposed beamforming update (7) and a full power transmission 
achieve the maximum rate improvement over the balanced solution 
(r*, r*) , i.e., reach the Pareto optimality with max-min fairness. 
This concludes the proof. 
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